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1. Introduction 

Stratospheric trace constituents may be detected through 
measurements of their infrared vibrational- rotational lines, and 
absolute column abundances can be determined from the measured line 
equivalent widths. If the individual line shapes are also measured, 
an inverse solution of the radiative transfer equation may be found 
which yields the density- altitude profile (for transient species) or 
the temperature- altitude profile (for a uniformly mixed gas). Measure- 
ments of individual line intensity profiles are very difficult by 
conventional spectroscopic methods, particularly when the gas is 
distributed throughout atmospheric regions having widely different 
total pressures. The line profiles will then consist of pressure 
broadened wings and doppler cores, imposing severe restrictions on the 
spectral resolution needed to fully resolve the line. Infrared 
heterodyne spectroscopy (iRIiS) is well suited to atmospheric studies 
because of its high specific detection sensitivity and its ultrahigh spectral 

resolution (~ 5 MHz = .00017 cm ' l ). At this resolution, individual 
line shapes can be measured without significant instrumental distortion 
and the analysis may be carried out in a straightforward manner. 

Application of heterodyne spectroscopy to stratospheric studies 
is still in its developing stages, and only preliminary attempts at 
measurements and some feasibility studies have been made (Seals and Peyton, 
1975; Seals, 1974; Peyton et al. 1977; Menzies and Seals, 1977). The 
purpose of this paper is to present additional calculations on the 
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feasibility of IRHS stratospheric composition measurements, and to 
report preliminary measurements and analyses of C0 n and 0 absorption 
lines measured with ground-based IRHS systems. 

Stratospheric measurements may be carried out either in absorption 
(observing the sun or the moon), or in self-emission (observing 
atmospheric emission lines directly). Absorption measurements are 
intrinsically absolute and the precision is limited by the signal- to- 
noise ratio (SNR) achieved. For solar measurements, the SNV may be 
as high as ~100,000, depending upon the volume mixing ratio of the gas, 
the strength of the line, the integration time, and the spectral bandwidth 
chosen. Lines with optical depths as small as 10 can be measured in this 
mode. The expected SNR's for self-emission are much lower, thus species 
investigated in this mode are limited to those having nearly optically 
thick infrared lines. 

This paper is limited in scope to ground-based solar- absorption 
measurements. The basic theory employed in making the radiative transfer 
calculations is discussed in section 2. Section 3 deals with the feasibility 
and application of an IRHS for stratospheric measurements, treating 
the detection sensitivity, minimum detectable volume mixing ratios of 
stratospheric constituents, and line shapes of infrared absorption lines 
in the stratosphere. The method of infrared spectral line inversion 
is briefly discussed in Section 4. Re^'ilts of solar-absorption 
measurements and analyses of C0 o and 0^ lines made with COp - laser 
and tunable-diode- laser heterodyne spectrometers are given in Section 5. 
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described by Selby and McClatchey (1972), Snider (1975) and Goldman and 
Snider (1975). 


The element of column density (du) may also be written in terms 
of the mass mixing ratio 

du = -^ q^(P) sec 0 dP, (5) 

where q^(P) - qY(P) rai /m^, with the masses referring to the species (i) and to 
the average molecular mass of the gas. The absorption coefficient at 
frequency v for species i is obtained by summing the contributions from 
all lines 

k = E k (P,T) , (6) 

where the absorption coefficient for an individual line may be written 


k^(P,T) = b dv = S t (T) b(P,T) , 


(7) 


where S^ is the integrated line intensity in cm ^ (cm atm) ^ and b 
is the line shape function. 

A convenient technique for calculating the absorption coefficients 
for an individual line is based on separating the contribution from each 
line into a wing contribution and a direct contribution near the line 
center (v ^ v c + 6), where 6 is model dependent and in the present case 
is chosen to be cm ^ (Kunde and Maguire, 197*+). The wing contribution 
is calculated by assuming a Lorentzian line shape, and the direct contribu- 
tion by using the mixed Doppler- Lorentz line shape (Voigt line shape) 
formula for which the function b is 
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the observed atmospheric spectrum (0.25 cm ^ resolution) is given 
in Fig. 3 (Kunde et al. 1977). 

(a) Sensitivity of a Heterodyne Spectrometer 

In a heterodyne receiver, the infrared signal from the source 
is mixed with coherent radiation from a local oscillator (such as a 
C0 2 laser or a tunable diode laser) and the difference frequency signal 
is detected by radio frequency techniques. The basic theory and 
sensitivity considerations of heterodyne spectroscopy have been 
given in the literature, (e.g., Teich, 1971; Abbas, et al. 1976). 

Here, we consider its sensitivity for detection of stratospheric 
molecules . 

The SNR of a heterodyne receiver is 

1/2 


S 

N 
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where A is the total system degradation from the ideal quantum detection 

limit, t is the integration time, B is the bandwidth of a single 

resolving element of the system, and P is the source power contained 

s 
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in the bandwidth and the etendue (AfJ ~ \ ) of a heterodyne receiver. 
Equation (13) yields an expression for the noise- equivalent radiance 
of a heterodyne receiver as : 

(HER). 
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Equation (14) indicates that the sensitivity of a heterodyne receiver 
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improves rapidly with longer wavelength, varying as \ . In principle, 

the spectral bandwidth (B) may be decreased arbitrarily for increasingly 
higher resolution, but only at the expense of a corresponding decrease 
in sensitivity. 

We now estimate the minimum detectable column density and volume 

mixing ratio. We assume a homogeneous isothermal atmosphere at 

v 

temperature T with uniform vO ime mixing ratio q^ and source surface 
temperature T . Equations (1) and (2) then give an expression for 

the minimum detectable vertical column density for the solar absorption 

mode 


(NER) 


Het 


U min k (T,P) B (T ) sec 9 * 
v v ® 


( 15 ) 


(16) 


The minimum detectable volume mixing ratio then becomes 

(NEd) , 

v _ x 7 Het 

Tnin k^(T,P) B^(T^) u^, sec 0 

where u^ is the total vertical column density and is the solar 
brightness temperature. 

For a moderately strong atmospheric line (S ~ 1 cm 1 (cm-atm~), 
k ~ 50 (cm- atm) 1 and for instrumental parameters of B = 5 I 'Hz, 
x = 1000 sec, A - 30, sec 6 = 1, the noise equivalent radiance and minimum 
detectable volume mixing ratio are 


(NER) 


Het 


and 


8 x lo" lS ( W cm" 2 Sr" 1 Hz" 1 ) 
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q v . ~ 2 x 10 or 2 pptv 
man ** 


( 17 ) 

(18) 


This detection limit has been derived only for 
obtaining an order of magnitude estimate of the sensitivity of an IRHS. 
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^^6* indic&t6s that at tills sensitivity almost till of* the important 
stratospheric constituents may be detected with a 
heterodyne spectrometer. Note that while Eq. (18) is based on a 
path length of one air-mass, slant paths at balloon altitudes can 
increase this to ~ 20 air-masses with a concomitant increase in 
sensitivity . 


(b^ Measurements of Stratospheric Absorption Lines 

The calculated transmittances for the mid- latitude summer 
model a tmosphere oi Selby and McClatehey (1^7*2) and the volume mixing 
ratios of Fig. 2 are given In Figs. 3-o. The computed spectra exhibit 
absorption lines of C0 o ,0 , OH, and N.^0. The water vapor concentration 
of the model atmosphere was modified to yield a total vertical column 
density of o8f n-atm, corresponding to field conditions during the 
observations discussed in section 5. 

An important consideration in the measurement- of line profiles and 
in the subsequent, retrieval of information is the tyre of lineshape 
involved. 1’he widths are on the order of 1 GHz with narrow cores. 

The widths of moderately strong lines of primarily stratospheric 
molecules (e.g. 0 . ,LC0) are on the order of a few hundred MHz. 

Representative lineshapes have been computed for stratospheric (z = 4s Km'' 
and tropospheric (z = 12 km and z = 0 km) levels for the a4o.55 cm" 1 (p?4) 
line of the ^k>l cm 1 C0 2 band and are shown in Fig. 7. At 4*5 km the doppler 
half-width is ~ 8 x 10 cm (24 MHz) and the core of the line is 
essentially shaped by doppler broadening. For Av greater than 
^l.o x 10 ' cm (50 MHz) the lineshape is dominated by pressure 


1 .? and 0 km levels the entire line is shaped 


broadeni ng . At. the 

by pressure broadening effects. Measurer ents 'f deppler cores v formed 

in the stratosphere^ thus require high spectral resolution (~ *' m?.' wheroa: 

lower resoluti »n MHc ^ will he sufficient in th.e wings (formed in the 

troposphere' 1 . 

k. Inversion of Infrared dpec t. ml hines 

Atmospheric po.ramet ers may be extracted by analysis o! individual 
spectral lines mens-.u-ed at sufficiently high spectral resolution tlirough 
an inverse solution of the radiative transfer equation {1\ The 
inversion of an appropr 1 at e.'.y oh*, sen lmie an provide cons . t uent. 
concentration profiles n both the stratosphere and the iruvsphere. 

Atmospheric sol f- emission it'. Kqn. (- 1 N usually may be ignored 
for solar observations, and the observed intensity can be written 


1° 1 ' h I-,' \' V 

v v 1 v 

obx .. 

The observed transmittance (t then is. 


ob ob v 

v v o . •> \ 

T ' V • 


Combining F,qs. (.'> ■ -1 the atmospheric transmittance may be 


written in the fo r. 


U\ r* " Z, it ,4 (T,F) q 4 (r) dr. 

V ... • V- 


Equation v 1 ' suggests Uvat r.he mixing raid os of the absorbing gas 
may be iteratively evaluated, with ^ used as a weighting function. 
The iterative scheme used here is similar to that, described in the 
literature (‘deals, 1- *7U ; dmith. l‘»70; Men. lies and dhahine, 1 


The plots ,'f weighting functions (normalised absorption 
coefficients) for some absorption lines of and iP,0 are 

shown in Mg.. (8)- (10). The weighting functions for 0 } in the 
d.W band reach peaks upto heights of - 30-39 *» Tor ground-based 
observations . The weighting functions reach peaks upto hr. 
for Cl!| ( and N ,0 absorption lines in ttie l.'.'O om -o 1- 
range. The altitudes which »y be probed by this method are 
determined by the shapes and the levels at which the weighting ‘''••notions 
ik ) reach their peaks. At altitudes whom -rvo-.-e oaden.i .. 


dominant, the Voight profile reduces to a i .'rente ian and 


K (P) * 


sor' 


(v-v 0 V + 




Setting the derivative of ^(P) with respect to a(P,T) equal to ,ero 
and solving for a(l \ indicates that for transitions with small values 
of lower state energy E" vso that the dependence of line strength on 
T may be ignored), k v (P) is maximum at levels where - i "v 0 ). 

At altitudes where doppler broadening is predominant (a o5 ’ 
k increases monotonically and thus cannot bo used as a weighting •.unction. 

The Inversion technique discussed here is thus .u.c.ui 

altitudes at which the lines are pressure broadened. A different mode of 
observation and inversion technique has to be employed for probing 
higher regions. 

Inversion of a specific line is carried out in the following way. 

The atmosphere is divided into N layers (35 in the present work'. 

The vertical temperature profile l\s' is assumed or evaluated, by 
analyzing the absorption line of a well mixed gas such as CO, (•’(*' 
need not be known as accurately for solar absorption measurements as 
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for the self-emission mode). An initial guess for the mixing ratio 

profile ,m is made (j denotes the level and m the iteration). The 

absorption coefficients and transmittances k J,m and t " 1 for an appropriately 

v k *k 

chosen set of frequencies v,. are calculated monochromatically by summing 
the contributions of individual molecular absorption lines as discussed 
in section 2. An improved mixing ratio profile q<3 i s then estimated 
at the levels corresponding to the peaks of the weighting functions from 

<r ’ = qr* (£n t /in t ) (23) 

V V 

The complete profile is interpolated from the corrected values 

and the iterative process is continued until the calculated t™ converge 

to the observed values of t°‘' over the entire set of chosen frequencies 

v 

in a least squares sense. The process is stopped when the RMS 
residual is of the same order as the noise level of the measurements. 

The weighting functions are chosen such that a sufficient number 
of levels are included for interpolation of the gas mixing ratio 
profile. The vertical resol i ion of the retrieved profile is then 
dependent on the shapes of the weighting functions, which are determined 
by the line strength and the lower state energy of the transition. 

If the weighting functions are sharply peaked, the .nversion procedure 
converges in 15-20 iterations. 

The mixing ratio profiles retrieved by inverting synthetic (O^CH^, 

NgO) and measured (CO,,) absorption lines (Figs. 4-6,14) are shown in 
Figs. 11 and 12. The initial guess profiles and the model profiles 
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are also shown. These results demonstrate the internal consistency 
of the method, and indicate that inversion of infrared absorption lines 
measured with sufficiently liigh spectral resolution can lead to vertical 
concentration profiles of stratospheric constituents to heights ~ 30 km 
with vertical resolution ~ 4 km. 

5. Atmospheric Absorption Measurements 

The atmospheric absorption measurements reported here are based 
on solar observations of absorption lines of C0 o and 0^ in the 10p,m 
band. The measurements were made with a C0,,-lascr heterodyne spectro- 
meter at Goddard Space Flight Center, and a tunable diode- laser heterodyne 
spec treme lex at MIT. 

A schematic of the Goddard CO,-, -laser spectrometer is shown in 
Fig. 13. Tliis system (Mumma , et al. 1977) j presently inter- 
faced with a 48" telescope at the optical test site in Ureenbelt, MD., 
employs a CC^ laser as local oscillator, a HgCdTe photodiode as a 
mixer^and two sets of IF filter banks. The high resolution filter 
bank consists of 40 sequential filters of 5 MHz resolution covering 
a total of 200 MHz bandwidth. The low resolution filter bank has 
23 filters of 50 MHz resolution with a total bandwidth of 1 GHz and 
one filter of 1 GHz bandwidth. The C0 o laser is grating tunable 
from line to line, so that the detection of atmospheric lines 
is limited to those which lie within + 1.5 GHz of a CO^ laser line. 

This restriction requires knowledge of the line center frequencies 
of the transitions to be observed. 


12 


Fig. 14 shows the observed profile of an atmospheric CO^ line 
v o = 9^7.7072 cm 1 ) measured on Feb. 10, 1977 with the low resolution 
filter bank. The Doppler broadened core, formed at high altitudes, can 
be identified. This line was inverted and the retrieved mixing ratio 
profile for COp is shown in Fig. 12. The computed absorption line 
based on the retrieved mixing ratio profile is shown by the solid 
curve in Fig. l4. The retrieved profile shows that C0 y is well mixed to 
at least 50 km, in agreement with independe t results, (e.g. Reiter 
et al., 1975). 

Fig. 15 shows an 0^ absorption line (v Q - 1044.9787 cm near 
the C0 2 - P(22) line in the 19.4vim band) measured on April 8, 1977 
at various values of air-mass (sec p). The high and low resolution spectra 
are shown on the le'’t and right hand sides respectively. With increasing 
air-mass, decreasing intensities and increasing line saturation are 
evident in both spectra. These measurements are of a preliminary nature. 
This line, with lover state energy of 707.21 cm ^ and lire strength 
S - 0.011 cm ^ (cm- atm) is not well suited for inversion. However, 
a measurement of this abs^ option line profile as a function of air-mass, 
indicates the capability of 1RHF> for investigating the daytime variation 
of 0^ concentration in the stratosphere. The nightime variations may be 
measured by making similar observations of the moon. 

Atmospheric solar absorption measurements were made with a 
tunable diode- laser spectrometer in Cambridge, MA. , on Jan 2, 1977 
around noon with the solar zenith angle between 67° to 64°. The 
air temperature was ~ 272 °K with relative humidity ~ 50%. 
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A detailed description of the tunable diode laser spectrometer 
and preliminary results have been given elsewhere (Frerking and 
Muehlner, 1977). The spectrometer employed a liquid nitrogen cooled 
FbSnSe diode laser as a local oscillator and a HgCdTe photodiode as a 
mixer. The diode laser was continuously tunable over the spectral 
range of 1010. O to 1011.8 cm \ The total bandwidth (2B) of the IF 
filters used in these measurements could be chosen as 200 MHz or 70 MHz, 
providing spectral resolutions ot' 0.0067 cm 1 or 0.0023 cm 1 respectively. 

The atmospheric spectrum taken with a spectral resolution of 
0.0023 cm ^ on Jan. 2, 1977* shows absorption lines of 0, in the 
spectral range of 1010.9 to 1011.8 cm 1 (Fig. 16). The spectral 
data on the eight 0, lines identified in Fig. lob are given in Table 
II. Since the atmospheric transmittance was not directly calibrated, 
an absolute transmittance scale was established by normalizing to the 
value at 1011.2 cm ^ calculated from a mid- latitude model atmosphere 
with a total water vapor content 'f 625 cm-atm (preeipi table water = 0.50 
cm). This water content corresponds to local conditions at the time 
of observations (Radiosonde data. National Weather Service). 

The strong ozone lines with S ~ 0.3 (in cm-atm units) are saturated 
(Fig. lob). The saturation is particularly strong because the air- 
mass during the measurements was ~ 2-3 (solar zenith angle 65°) 
even though the measurements were made around noon. The saturation 
results in a loss of spectral information from the upper stratosphere, 
since the weighting functions for frequencies near line center peak 


at heights — 25-35 km. A line of moderate intensity has to be used to 
probe the upper stratosphere from the ground. In addition, transitions 
with smaller values of LC energy are needed to reduce the temperature 
dependence. A comparison of the numerical values of 3 and LS 
energy listed in Table III, shows that lines 3j4, and 8 are most 
suitable for analysis. Lines 4 and 8 are saturated, but 
line 3 with a lower value of 3, may be used in conjunction with line 
4 or 8 for analysis. The weaker .lines such as 1 oi 2, with larger LS 
values, have weighting functions which do not have sharp peaks and 
are thus not suitable for inversion. 

Lines 3 and 8 , with line centers at = 1011.1039 cm ^ and 

v Q = IOH .6636 cm ^ respectively, are shown at high resolution in 

Fig. 17a ,b. The retrieved mixing ratio profile of 0- obtained 

through an inversion of these two lines is shown in Fig. 18. Because 

of saturation of the two lines ar.d insufficient resolution near the 

line center, the retrieved profile is accurate up to heights of 

-6 

about 25-28 km only. A fixed value of q v = 2.0 x 10 is assumed 
at a height of 50 km and the complete vertical profile is then 
interpolated. The retrieved profile corresponds to a total vertical 
O^ column density ~j 0.37 cm atm. This value compares with a column 
density of 0.398 cm-atm of' the mid- latitude model 0 ^ profile shown 
in Fig. 18. 

The computed transmittances for the two lines based on the 
retrieved 0^ profile agree well with the observed values (Fig. 17). 

The synthetic spectrum was calculated over the spectral range 
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1019.9”10U-*8 cm 1 from the retrieved ozone profile (Fig. 16c). 

The line parameters such as line centers, line strengths and 
lower state energies were taken from the AFCRL line parameter 
atlas (McClatchey et al., 1973)* The half-width for ozone used 
was 0.07 cm ^ at NTP. No attempt has been made to correct the 
line centers in the synthetic spectrum to conform with the observed 
positions. 

A comparison between the observed and synthetic spectra (Figs. 

16b and loc) shows that there are significant differences, indicated 
by missing or extra lines. No attempt lias been made to 
achieve a closer match between the two plots, but the differences are 
very likely due to inaccuracies in line parameters used in generating 
the synthetic spectrum. Existing uncertainties in the line centers 
and intensities are evident when comparing y^ line parameters of 
Barbe et al ( 1977 ), Menzies ( 1976 ), and the AFCRL atlas (McClatchey 
et al., 1973 ). Extra lines could also be due to other gases. 
Conclusions : 

The measurements and analyses presented in this paper indicate 
the usefulness of heterodyne measurements in the infrared for detec cion 
of a large number of stratospheric constituents. The vertical mixing 
ratio profiles and the total abundances of some constituents may be 
determined with high vertical resolution (~ 4 km) and the diurnal 
•variation may be followed by tracking tlie sun or the moon. Interpreta- 
tion of heterodyne measurements and the accuracy of anulysis is limited 
by the availability of accurate line parameters. A parallel program 
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Table I. Some Selected Gases with Transitions 
In the Intermediate Infrared 


Gas 

Accessible 
Spectral Region (cm 


°3 

1000-1100 

* 

NO : 

NgO 

1190,1290 

* 

N0 £ 

1621 


HNO 

890 


V > 5 

1000 


HO. : 

HgO 

900-1200 


ho 2 

1101,1390 


%°2 

880,1266,1380 


Cl : 

CCl^ 

768,789 


CFC<t^ (Freon 11) 

847 


CF 2 C ^2 ( Freon - 12 ) 

920-930 

* 

CH^C-t 

732,1355 

* 

CIO 

850 

* 

C-tONOp 

780,1490 


Br: 

CH-.Br 

954,1443 


BrO 

853,1^83 


NH^ 

930-970 


ch 4 

1306 


CH 3 F 

1047 

* 


*Line data known accurately for at least some lines. 


(Hot Bands) 



Table II. Spectral data for the ozone ataorption lines in 
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